ImgList = {
    rgb2ycbcr(image)
    rgb2ycbcr(image_b)
    rgb2ycbcr(image_bh)
    rgb2ycbcr(image_h)
    };

ImgList = cellfun(@(x)x(:,:,1), ImgList, 'uniformoutput', false);

%%
r = ImgList{3} - ImgList{2};
y = ImgList{end} + r;% / sum(r(:));

close all
imshows(r)
imshows(ImgList{end})
imshows(y)

srimg.psnr(ImgList{1}, ImgList{end}, 3)
srimg.psnr(ImgList{1}, y, 3)
